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ABSTRACT 

We use a sample of 8298 galaxies observed as part of the HST ifi6o-band GOODS 
NICMOS Survey (GNS) to construct the galaxy stellar mass function both as a func- 
tion of redshift and stellar mass up to z = 3.5. Our mass functions are constructed 
within the redshift range z = 1 — 3.5 and consist of galaxies with stellar masses of 
M* = 10 12 M Q down to nearly dwarf galaxy masses of — 10 8 5 Mq in the low- 
est redshift bin. We discover that a significant fraction of all massive M* > lO n M0 
galaxies are in place up to the highest redshifts we probe, with a decreasing fraction of 
lower mass galaxies present at all redshifts. This is an example of 'galaxy mass down- 
sizing', and is the result of massive galaxies forming before lower mass ones, and not 
just simply ending their star formation earlier as in traditional downsizing scenarios, 
whose effect is seen at z < 1.5. By fitting Schechter functions to our mass functions 
we find that the faint end slope ranges from a = —1.36 to —1.73, which is significantly 
steeper than what is found in previous investigations of the mass function at high red- 
shift. We demonstrate that this steeper mass function better matches the stellar mass 
added due to star formation, thereby alleviating some of the mismatch between these 
two measures of the evolution of galaxy mass. We furthermore examine the stellar 
mass function divided into blue/red systems, as well as for star forming and non-star 
forming galaxies. We find a similar mass downsizing present for both blue/red and 
star-forming/non-star forming galaxies, and further find that red galaxies dominate 
at the high mass end of the mass function, but that the low mass galaxies are mostly 
all blue, and therefore blue galaxies are creating the steep mass functions observed 
at z > 2. We furthermore show that, although there is a downsizing such that high 
mass galaxies are nearer their z = values at high redshift, this turns over at masses 
M* ~ 1O 1O M0, such that the lowest mass galaxies are more common than galaxies at 
slight higher masses, creating a 'dip' in the observed galaxy mass function. We argue 
that the galaxy assembly process may be driven by different mechanisms at low and 
high masses, and that the efficiency of the galaxy formation process is lowest at masses 
M* ~ 10 10 Mq at 1 < z < 3. Finally, we calculate the integrated stellar mass density 
for the total, blue and red populations. We find the integrated stellar mass density of 
the total and blue galaxy population is consistent with being constant over z = 1 — 2, 
while the red population shows an increase in integrated stellar mass density over the 
same redshift range. 

Key words: galaxies: evolution-galaxies: general-galaxies: luminosity function, mass 
function 



1 INTRODUCTION 

A deep understanding of the high redshift universe is vi- 
tal in order to complete our knowledge of galaxy formation, 
and hence uncover the history of the universe as a whole. 



With the continued development of instrumentation and of 
new analysis techniques our knowledge of the high redshift 
universe is increasing rapidly. Thanks to new deep imag- 
ing and multiobject spectroscopy we now have the power 
to routinely look back at the universe over cosmic time 
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to witness step by step evolution. The result of this has 
been a wealth of observations of large samples of galax- 
ies over a large redshift range in vario us surveys, such as 
within GO ODS iGiavalisco et ai1l2004l), COSMOS a nd z- 
COSMO S dScoville et alj|2007al and iLillv et alj 12007ft and 
AEGIS ijDavis et alJl2007t ). These large samples give us the 
power to achieve statistically meaningful results concerning 
the evolution of galaxy properties. 

As a result a detailed picture is beginning to 
form regarding when galaxy stell ar mass is bu i lt up 



over cosmic time . 
Drorv et all 120051. 



Various studies IIDickinson et al 
IConselice et al.ll2007l, lElsner et al.1 



2003 



2008 



Perez-Gonzalez et al.l I2008T ) have focused on the change of 



stellar mass density with time, and have seen generally con- 
sistent results. These studies find that the integrated stellar 
mass density decreases at higher redshift, as expected since 
the ongoing process of star formation increases the amount 
of stellar mass in the universe over time. These studies also 
show that roughly 50% of the mass density of the universe 
is in place by z ~ 1. This implies prior to this redshift the 
comoving stellar mass density has a more rapid evolution. 
This ties in with studies that show that the star formation 
rate peak is at Q 1 - 2 (|Madau et al.lll99fj , iHopkindliooi. 
iHopkins fc Beacomll2006T ). 

While studies of the star formation history has been the 
traditional way to examine and probe galaxy evolution, a 
great deal of research has been carried out investigating the 
evolution of galaxies using the galaxy stellar mass function. 
Early examples of the measurement of the local g alaxy mass 
functi on at z ~ have bee n carried out by e.g. I Cole et ail 
i|200ll ') and lBell et all l|2003h , using surveys such as 2MASS, 
2dF, and the Sloan Digital Sky Survey, data. These and 
other investigations construct the stellar mass function, and 
integrated stellar mass density within the local universe, and 
hence provide a vital benchmark for comparison with higher 
redshifts. 

It is important to understand the stellar mass func- 
tion of the local universe in itself as this traces the inte- 
grated star formation and mass assembly history over the 
entire universe. However, by extending similar studies to 
higher redshift we can investigate not only galaxy growth 
as a function of stellar mass, but also the growth of inte- 
grated galaxy stellar mass with time. That is, we can trace 
the evolution of stellar mass for galaxies of different stellar 
masses, i.e., very massive galaxies vs. lower mass galaxies, 
over time. This gives us insights into the growth of galax- 
ies of different stellar masses due to star formation, mergers 
and other assembly processes. Many studies have investi- 



large area surveys dFontana et al. 2004, Bundv et al. 


2006, 


Borch et al.ll2006l. iFranceschini et al.1 120061 iBell et al. 


2007. 


Bolzonella et al. 


2009). This has been further extended to 



even higher redshifts u sing deeper surveys, g e nerally within 
a much smaller area |Conselice et al.l 120051. iFqntana et all 
2006), iPerez-Gonzalez et alj|2008t lKaiisawall2009l ). 



These studies of the stellar mass function at high red- 
shifts allows us to form a picture of the high mass end 
of the stellar mass function as a function of redshift. It is 
widely agreed that the dominance of star formation within 
massive galaxies ends much e arlier than within low stel- 
lar mass galaxies by z ~ 1 (|Bundv et al.l l2006h . This is 
one form of galaxy downsizing, whereby the higher stel- 



lar mass galaxies have their star formation truncated, or 
gas depleted, earlier than lower mas s galaxies. This for - 
mation scenario was first observed by ICowie et al. I (jl996l ) 
an d subsequently obs e rved in various stu d ies including 
e.g.lBauer et al. I i2005h.lFeulner et al. I (|2005l ). lBundv et all 
(2006), and IVergani et al. I (|200sh . Downsizing is now ac- 
cepted as part of the formation scenario of galaxies, but is 
not yet fully understood. What is also not fully understood 
is whether the related, but different process of mass downsiz- 
ing, is occurring in tandem, such that the high mass galaxies 
form their stellar mass earlier than lower mass galaxies, and 
when this stellar mass differentiated galaxy formation pro- 
cess first reveals itself. 

Despite the considerable work done investigating the 
high mass end of the stellar mass function, there are still 
many issues that are not yet fully understood. Firstly, the 
generally shape of the high redshift stellar mass function is 
not well described. Nearby stellar mass functions are well 
fit by a form of the Schechter function, and there are vari- 
ous investigations regarding how the parameters of this fit 
change with redshift. Secondly, difficulties with obtaining 
deep data at high redshifts mean that the low mass end of 
the stellar mass function has not been explored as fully, or to 
as high a redshift. More recent work has started to uncover 
a possible steepening with redshift of the low stellar mass 
end, and a "dip " in the intermediate stellar mass range (e.g. 
iKaiisawal 12009,) . It has been suggested that this is a result 
of evo lution of different galaxy populations dr i ven by their 
mass dBolzonella et al"]|2009l . iDrorv et al"]|2009l . Illbert et af] 
l2009bl . iPozzetti et alj|2009h . The exact nature and reasons 
behind such features of the stellar mass function are not 
well understood, and this can only be improved upon with 
deeper, more robust data. This will then lead to a better un- 
derstanding of how the populations of galaxies, as defined 
by stellar mass, change and evolve over time. 

In this paper we use data from the GOODS NICMOS 
survey (GNS) to investigate how the stellar mass function 
evolves from z = 1 to z ~ 3.5. By examining the stel- 
lar mass functions of galaxies ranging in stellar mass from 
M* = 1O 12 M to as low as M, = 1O 8 ' 5 M we investigate 
when, and which galaxies are forming at various epochs in 
the universe. The depth of the GNS data is such that we are 
able to probe over a factor of 10 3 the stellar mass evolution 
up to z ~ 3.5 and importantly trace how galaxies of differ- 
ent masses are evolving with time. We find throughout this 
paper a differential in the stellar mass function and how it 
evolves, revealing a strong stellar mass dependence in the 
galaxy formation process. We describe this and give some 
general explanations for how this differential evolution can 
occur due to different physical processes. 

The paper is set out as follows: Section discusses the 
GOODS-NICMOS Survey, the galaxy sample and how the 
data used in this paper was obtained. Section \3 . 1 1 examines 
the galaxy stellar mass functions of all the galaxies in various 
redshift bins. In Sections 13.21 and 13.31 we split the galaxies 
into blue and red and star forming and non-star forming 
respectively. Section[2]describes the calculation of the stellar 
mass densities for the total sample and for the red and blue 
populations. Sections [5] and [6] contain the discussion and 
summary of our findings respectively. Throughout this paper 
we assume !1m = 0.3, Q.a = 0.7 and Ho = 70 km s _1 Mpc -1 . 
AB magnitudes and a Salpeter IMF are used throughout. 
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2 DATA REDUCTION AND THE GNS 

2.1 The GOODS-NICMOS Survey 

The galaxy sample used in this work is taken from and 
imaged as part of th e GOODS-NICMOS Survey (GNS) 
|Conselice et alj|201ll ). The GNS consists of 60 pointings of 
the HST NICMOS-3 camera utilising a total of 180 orbits. 
The field of view of NICMOS-3 is 51.2 arcsec x 51.2 arcsec 
with a pixel scale of 0.203 arcsec/pixel. Each NIC3 tile was 
observed in six exposures, which combined give a pixel scale 
of 0.1 arcsec/pixel with a point spread function (PSF) of ~ 
0.3 arcsec full width half maximum (FWHM). Within the 
GNS we find 8298 galaxies in the H wo band (F160W), and 
the 60 pointings are designed to contain as many massive 
galaxies (M» > 10 11 Mq) as possible. These massive galaxies 
are in the redshift range z = 1.7 — 2.9 and at a depth of three 
orbits. The galaxies are redshift sel ected based on t h eir op - 
tical to infrared colour described in IConselice et all l|201lh . 
Details on the data reduction pipe line are described in 
iMagee. Bouwens fc Illingworthl (|2007h . The detections and 
photo metry were done using SExtractor l|Bertin fc Arnoutsl 
1996). At 5a the limiting magnitude is _f/i6o=26.8, a marked 
improvement on the GOODS ground based data at K =24.5 
|Retzlaff et al J 120091 ) . Further description of the GNS, the 
pointi ngs and the target selection is given in lConselice et all 
J201ll). Other ana l ysis o f the GNS da t a set can be found 



Buitrago et aJ~l j|200St). iBluck et all |2009l), Bauer et al 



(2010, submitted), IBluck et all (|2010h . iGriitzbauch et all 
|201ll ). 



2.2 Photometric Redshifts 

Thanks to the large amount of optical data covering the 
GOODS fields, the -ffieo band sources are matched to a 
catalogue of B,V,i and z band data. This photometry is 
available down to a limiting magnitude of B~28.2, and the 
matching is done within a radius of 2". However the mean 
separation between the optical and -H160 band coordinates 
is markedly better than this with (r) ~0.28±0.4", which 
is roughly the resolution of NICMOS. With this multiband 
data template spectra we fit to the BVizH photometric data 
points. This was handled in two different ways to overcome 
the degeneracy in colour-reds hift space. 

The first uses HYPERZ jBolzonella et all 120091 ) which 
is the standard \ 2 minimisation technique. HYPERZ uses 
model spec tra which were cons t ructe d using the evolution- 
ary codes of lBruzual fc Chariot! |l993t ) . Here we use five evo- 
lutionary types, E, Sa, Sc, Im and a sing l e bur st scenario. 
The reddening law is that of lCalzetti et al.l {2000). The most 
likely redshift is then computed in the age, metalicity, red- 
dening parameter space, giving the best fit redshift, corre- 
sponding probability and several other best fit parameteres. 

The second method used to obtain our photometric red- 
shifts is the Bayesian approach using BPZ. This is a similar 
fitting method to HYPERZ but employs empirical SEDs 
rather than model ones. As well as this, BPZ uses the max- 
imum likelihood in the same parameter space as HYPERZ, 
but with the addition of empirical information regarding the 
likelihood of a certain combination of parameters. This is 
known as prior information, or priors. Here, we used the dis- 
tribution of magnitudes of different galaxy t ypes as a func 
tion of redshift as the priors (from HDF-N 



Therefore in this case, not only does the code find the best 
fit redshift solution and spectral type, it also consideres how 
likely it is to find a galaxy of that type and magnitude at 
that redshift. 

A total of 906 spectroscopic redshifts are avail- 
able for compari son with our photometr i c reds hifts from 
the GOO DS-N jBarger. Cowie. fc Wand 120081 ) and the 
GOODS-S |Wuvts et alj|2008l ) fields. The reliability of pho- 
tometric redshifts is defined as Az/(1 + z) = (z spec — 
z p hot)/(l + z spec ), the median error ({Az/(1 + z)}) and the 
r.m.s (a) are as follows. For HYPERZ, the GOODS-N gives 
(Az/(l + z)) = 0.027 with cr=0.04, and for the GOODS-S 
(Az/(l + z)) = 0.043 with cr=0.04. For BPZ, the GOODS- 
N (Az/(1 + z)) = 0.07 with cr=0.05 and for the GOODS-S 
(Az/(l-\-z)) — 0.07 with cr=0.06. In both codes the extreme 
catastrophic outliers, as defined by |Az/(l + z)\ > 0.5, are 
around ~6%. The dependence of the redshift on the -H160 
magnitude shows the reliability of redshift within our selec- 
tion methods and in this case HYPERZ sh ows a slightly bet- 
ter performance ijGriitzbauch et al]|201l[ ). Our work shows 
good agreement with photometric redshifts from past sur- 
veys (e.g. FIREWORKS), although our sample is bright, 
and it is unclear if this would be the case for fainter galax- 
ies. 

We further investigate the performance of HYPERZ at 
different redshifts, at low redshift (z < 1.5) and in the red- 
shift range of 1.5 ^ z ^ 3, which is the redshift range of 
the galaxy sample we use in this study. For the high redshift 
sample we obtain an average offset (Az/(1 + z)) = 0.06 and 
a RMS of (7az/(i+z) = 0.10, with a fraction of catastrophic 
outliers of 20%. Here catastrophic outliers are defined as 
galaxies with |Az/(l + z)\ > 0.3, which corresponds to ~ 
3 times the RMS scatter. Galaxies below z = 1.5 show a 
slightly lower, but still comparable scatter of (7a z /(i+z) = 
0.08, however the outlier fraction decreases dramatically to 
only ~ 2%. Furthermore, we simulate the effects of the pho- 
tometric redshift errors on our results throughout this paper. 



2.3 Stellar Masses and Colours 

The stellar masses we use are obtained from the same 
BVizH catalogue as used to measure the photometric red- 
shifts described previously, using a stand ard multicolour 
stella r po pulation fitting tech nique (e.g., IConselice et al.l 



120081 and IConselice et alj [201ll 'l. The photometry is fit to 
different star formation histories, based on the redshift of 
the galaxy, with spectroscopic redshifts used when available. 
This produces a distribution of likely stellar masses, rest 
frame optical colours and various other parameters based 
on a Baysian approach. The model SEDs are constructed 
from lBruzual fc Chariot! |2003l ) models, and the star forma- 
tion history is characterised by an exponentially declining 
model with various ages, metalicities and dust extinctions. 
The star formation rate is parameterised by an e-folding 
time and an age such that, 



SFR(t) ~ SFRo xe~~. 



(1) 



t ypes as a tunc - 
■ iBem'tezl feoOOy 



The parameters in Equation [T] are varied randomly within 
the ranges; r=0.01 to 10 Gyrs, t=0 to 10 Gyrs, metallic- 
ity= 0.0001 to 0.05 and the dust content is parameterised 
by r„=0.0, 0.5, 1, 2. These model SEDs are then fit to the 
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observed photometric data points using a Bayesian approach 
resulting in a likelihood distribution of stellar mass, age and 
absolute magnitude for each possible star formation history. 
The stellar mass is determined based on this distribution, 
where the most likely stellar mass produces a peak in the 
distribution, and the uncertainty is the width. The final er- 
rors produced are a result of the models used and are found 
to be in the range 0.2 to 0.3 dex. While several of the other 
parameters produced by this method ar e not reliable, the 
stel lar masses and colours are robust fsee lBundv et al1l2006l 
and lConselice et al.l l2007 for further explanation.). It is also 
possible that the stellar masses are an over estimate due to 
the poor treatment of the TP-AGB phase in a star's life. 
The effects of this are much less important at the rest frame 
wavelengths used in this study, especially in the infrared H- 
band. By using the newer models of Bruzual and Chariot 
(2010) which include a more proper treatment of TP-AGB 
stars we find it only lowers the masses in our massive galaxy 
sample by <0.07 dex. This effect is much smaller than the 
stellar mass error and the effects of cosmic variance and it 
is therefore considered negligible in this work. 



2.4 Star Formation Rates 

To obtain star formation rates (SFRs) for our GNS galaxies 
we use the rest-frame ultraviolet (UV) luminosity, which we 
derive from observed rest frame optical light following the 
procedure described in Bauer et al. (2010, submitted). The 
UV luminosity is closely related to the level of ongoing star 
formation because it is mainly produced by short-lived O 
and B stars. An advantage of using UV light to estimate 
star formation is that it re mains largely u naffected by the 
age-metallicity degeneracy IjWorthevI I1994T ). but as is well 
known, star formation rate estimates derived from UV light 
are strongly affected by dust extinction. 

The 2800A rest-frame luminosity is calculated from the 
observed optical Hubble/ ACS z-band flux density, which 
corresponds to rest-frame wavelengths of 3400 - 2125A for 
z = 1.5 — 3 galaxies. To derive SFR^v we apply a sim- 
ple K-corre ction derived from the redshift of e ach object as 
log(l + z) (iKim et al.lll996l . iDaddi et al]|2004h and use the 
iKennicuttl jl998h conversion from 2800A luminosity to star 
formation rate assuming a Salpeter IMF: 



SFRuv (M s yr _1 ) = 1.4 x KT 28 L v (ergs s -1 Hz" 1 ). 



(2) 



Before dust extinction is taken into account fSection l2.5[) . we 
find at z = 1.5 a limiting SFRc/y = 0.28±0.1 Moyr" 1 and at 
z = 3.0, we find a limit of SFR UV = 0.98±0.3 M s yr _1 . The 
errors take into account photometric errors and the error in 
the conversion from a luminosity to a star formation rate 
(see also Bauer et al. (2011, in prep.). 



near-ultraviolet band and z ~ 2 galaxies show that the UV 
slope from the local starburst relation can be used to re- 
cover the dust attenuation o f a vast majority of moderately 
luminous galaxies at z ~ 2 |Buat et al. 120051 . ISeibert et al.l 
120051 . iReddv et alj|2010h . 

To determine the UV slope, we use the SED-fitting 
procedure. We fit a spectral energy distribution (SED) to 
multi-wavelength observations from o ptical to mid-infrared, 
follow ing the procedure described in iPerez-Gonzalez et al.l 
{200$) and Bauer et al. (2011). Briefly, the UV through 
MIR SEDs obtained for all sources in the GOODS fields 
were fitted with stellar population synthesis models. Then, 
the best-fitting templates were used to get synthetic estima- 
tions of the UV emission at 1600A and 2800A. From the syn- 
thetic, model-derived UV luminosities at 160 A and 2800 A, 
we cal culate the spectral slope, /3. We use the lCalzetti et all 
(2000) law to derive A2800 from the UV spectral slope, 
which we then apply to the UV-derived star formation rates. 
Using this method we find an average extinction value of 
A2800 = 1-6 ± 1.2 mag. We find that all galaxies in the 
GNS exhibit a range in SFKuv,corr between 0.2 Moyr -1 
and ~2000 Moyr -1 . A full description of our star formation 
measure is provided in Bauer et al. (2011). 



3 THE STELLAR MASS FUNCTIONS 
3.1 Total Stellar Mass Function 

To construct galaxy stellar mass functions our data has been 
split into five redshift bins between z = 1 and z = 3.5. These 
are chosen to have roughly the same comoving volume to re- 
duce fluctuations in the total number of galaxies. The high- 
est redshift bin, z = 3 to 3.5, contains around 300 galaxies, 
thus we maintain a good sample size throughout. However, 
we exclude the final redshift bin from much of our analysis 
as, in that redshift range, the H-band is no longer sampling 
primarily optical lightand no longer provides a good mea- 
sure of the balmer break, thus we find larger errors in the 
measurement of stellar masses. 

We calculate the number densities, i.e. the number of 
galaxies per co-moving volume per mass interval, and plot 
the stellar mass functions for each redshift bin in Figure [T] 
Plotted for comparison in Figure [1] as a solid b l ack cu rve, is 
the local stellar mass function from lCole et alj l|200lf ). Also 
plotted is the vertical line that indicates the mass limit of the 
survey in each redshift bin (explanation of the calculation for 
this can be found in Section ^. 1.1|) . These lower limits are at 
a very low stellar mass c ompared with other work su ch as 
Marchesini et alj i2009h . IPerez-Gonzalez et all (2008) and 
Font ana et a ll l|2006h . who do not probe stellar masses as 
low as we do (see Section [5] for discussion of this). 



2.5 Dust Correction 

To obtain reliable star formation rates from the rest-frame 
ultraviolet, we need to account for the amou nt of light ob- 
scured by dust, which is a non-trivial problem. iMeurer et al.l 
(1999) demonstrated a correlation between dust attenua- 
tion and rest-frame UV slope, /3, for a sample of nearby 
starburst galaxies (where Fx ~ A^). Updated studies of lo- 
cal galaxies using the Galaxy Evolution Explorer (GALEX) 



3.1.1 High Mass Bias and Completeness 

The GNS is specifically intended to maximise the popula- 
tion of high mass galaxies, and hence the galaxy sample will 
contain more galaxies with M* > 10 1 Mq than expected for 
a randomly positioned survey. The high mass end of the stel- 
lar mass function is corrected for this by computing the ratio 
of high mass galaxies in the GNS pointings to the number of 
high mass galaxies in the total GOODS fields. By comparing 
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Figure 1. The galaxy stellar mass function between redshifts 2 = 1 — 3.5 in the GNS. The dashed green curve is the Schechter fit to 
the data. The black points are the data that has been fitted, and the red points are the data that has been left out of the fit due to 
incompleteness. The final rcdshift bin is represented by open circles as it is not included in the bulk of our analysis as discussed in Section 
13.11 The solid pink vertical lines show the theoretical mass limits of the GNS survey (see Secti on l3. 1.111 . Also included for comparion is 
the solid black line which represents the local galaxy stellar mass function of Col e et al. I fcOOlT) . The data used to construct these mass 
functions can be found in Table [T] 




Figure 2. The parameter results of the Schechter fits. Left: The results of the fitting for M* with all parameters free, Middle: The result 
of the fitting for a with M* held constant (see Tabled, and Right: The results of the fitting f or <§>* with M* held constant. For each 
panel the black circles are the results from this work, the red triangles arc Perez-Gonzalez ct al . (2008) results, the g reen squares are the 
iKaiisawal j2009l ) results (from three different SED models) and the pink stars (pink line for middle panel) are the lElsner et all (2008) 
results. For each panel the final redshift point is plotted as an open circle as it is not considered in the analysis. The parameters of the 
final redshift bin are represented by open circles as they are not included in the bulk of our analysis as discussed in Section 13.11 
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Redshift Range 


M* 


0(xlO" 4 ) 


Q 


1.0 to 1.5 


11.43 


6.01 ± 1.05 


-1.36 ± 0.05 


1.5 to 2.0 


11.43 


7.53 ± 1.23 


-1.19 ± 0.06 


2.0 to 2.5 


11.43 


3.52 ± 0.89 


-1.50 ± 0.08 


2.5 to 3.0 


11.43 


1.11 ± 0.36 


-1.89 ± 0.11 


3.0 to 3.5 


11.43 


0.89 ± 0.22 


-1.73 ± 0.09 



Table 2. The values of the parameters from the Schcchtcr fit. M* is the mean value from fitting with all parameters free. <j>* and a are 
the result of the Schcchtcr fit with M* held constant. 



this to the fraction of the area that the pointings covered, 
we obtain a correction factor for the overdensity of galaxies 
with M* > 10 n M©, finding a value of 3.05. The high stellar 
mass galaxies (M* > 10 11 Mq) in the redshift bins 2.0 to 
2.5 and 2.5 to 3.0 we re thus divided by 3.05 before fitting 
ijConselice et al.ll201ll ). 

Another problem in constructing the stellar mass func- 
tions is incompleteness at the low mass end. As mentioned 
previously, Figure [T] shows the mass limit of the survey. This 
is a purely theoretical mass limit calculated from the central 
wavelength of the NICMOS camera. We use this to calculate 
the rest frame band observed by NICMOS for each redshift 
bin. From this we calculate the mass to light ratios of a max- 
imally old stellar population, again for each redshift range. 
We combined this with the luminosity limit calculated from 
the limiting magnitude of the survey and the luminosity dis- 
tance for each redshift bin. This gave us the mass limit of 
the survey represented by the pink vertical line in Figure [1] 

Figure [1] shows that galaxies appear to drop out before 
the calculated mass limit. This is due to the pipeline de- 
tection being less sen sitive to galaxies at the low mass end 
ijConselice et al.ll201lf ). thus we do not include these points 
in our fits. To determine when this occurs we look at the 
residuals between the local mass function and the calculated 
number densities. A change in trend of the residual repre- 
sents the loss of completeness and we do not fit below this 
part of the stellar mass function. The points not included in 
the fit are the red points in Figure [T] 

3.1.2 Fitting the Schechter Function 

We fit our stellar mass functions, within the errors on 
the n umber densities, with a Schechter function (jSchechterl 
1 19761 ) of the form 

ct>(M) = 0* • log(10) • [i < M - M *>]< 1+t *> ■ ffl Ep[-10 (M - M * ) ]. (3) 

In this equation the parameter M* is the characteristic mass 
at which the stellar mass function turns over, a parame- 
terises the slope of the faint end of the stellar mass func- 
tion, and (f>* is the scale factor. The first fitting of the stellar 
mass functions we performed was done by leaving all of the 
three parameters free. For the redshift bins z = 2 — 2.5 
and z = 2.5 — 3 we find that the highest mass points 
(M* > 10 11 Mq) have extremely large errors. To help con- 
strain the high mass end we combine the highmass points 
into one point in these redshift bins. 

The left hand panel of Figure [2] shows the results of 
M* for the first fitting, where we find there is very little 
change in M* over redshift (except at z = 3 — 3.5). This 
is in good agreement with various other studies such as 



iPerez-Gonzalez et all (I2008T ) and lElsner et al.1 (120081 ) (these 
points can also be seen in the left hand panel of Figure [2}. 
We then repeat the fitting holding M* constant at its mean 
value (see Table [2} over the whole redshift range, which al- 
lows for a better constraint on a and <f>* . Some previous work 
hav e found a constant a up to z ~ 2 (e.g. lFontana et al 1l2004l 
and lBorch et al.ll2006l ). In this work we find no evidence of 
this, hence we only use a constant M* for investigating the 
evolution of a and <f>* at 1 < 2 < 3. The values of the fitting 
parameters are shown in Tableland also compared to pre- 
vious work in Figure [2] For this repeated fitting we do not 
combine the high mass points as before, as the high mass 
end of the stellar mass function is already well constrained 
by M*. 



3.1.3 Errors and Simulations 

We calculate errors on our number densities using two Monte 
Carlo simulations taking into account the lcr Gaussian mea- 
sured error on the stellar masses and accounting for errors 
on the redshifts. We first use the measured error on each 
stellar mass and compute a Gaussian distribution of simu- 
lated stellar masses, for each galaxy, between ±3<r of their 
measured error. Then a new stellar mass was randomly se- 
lected from this Gaussian distribution, so that we obtain 
a simulated stellar mass for each galaxy. We then recalcu- 
late the number densities in the same way as we did for the 
original catalogue of stellar masses, and the error due to the 
catalogue mass error, is the difference between the simulated 
number density and the original catalogue number density. 

For the photometric redshift errors we construct a cata- 
logue of simulated redshifts using the same Gaussian method 
we use to obtain the simulated stellar masses described pre- 
viously. We assume Gaussian errors on the photometric red- 
shifts, and do not consider the catastrophic outliers. We then 
calculate the luminosity distance for both the original and 
the simulated redshifts of each galaxy. We calculate the ratio 
of these two luminosity distances and assume this is equal 
to the ratio of stellar masses. By multiplying the original 
galaxy stellar mass by the luminosity distance squared ra- 
tio, we obtain a new simulated galaxy stellar mass. As before 
we then recalculated the number densities using this cata- 
logue of simulated stellar masses, then subtract the original 
stellar mass from the simulated stellar mass to obtain the 
error. We then add the stellar mass and redshift errors in 
quadrature, as well as the error due to Poisson statistics, to 
obtain the total error on each number density. These are the 
errors shown on Figure [1] 

We also simulate independently the effects of the mea- 
sured errors on the stellar mass function by randomly alter- 
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Table 1. 



The total galaxy stellar mass functions seen in Figure 



ing each stellar mass by ±0.25 dex. We chose this number 
to investigate the extremes of the effect of large errors. This 
shows how a mass bin containing less galaxies will be effected 
by the errors on the stellar mass more than a fuller mass bin. 
This is due to the Eddington bias, i.e. if a mass bin is rela- 
tively empty, galaxies from nearby fuller bins will spill over 
in to these bins more readily due to measurement errors. 
We reanalyse the altered masses and find that the largest 
variation is in the emptier, high mass bins. These variations 
lie within the calculated errors on the number densities and 
thus we can still reconstruct the same Schechter function 
parameters. 



3.1.4 Inspection of the Mass Functions 

There are several features that can be noted through in- 
spection of the stellar mass functions, and from the best fit 
parameters. Firstly, the massive galaxies (M* > lO n M ) 
are present, with a similar number density as at z — 0, up 
to a redshift z = 3. The low mass galaxies do not reach 
the local value of the number density until after the massive 
galaxies. This is a downsizing in terms of stellar mass over 
a large range. We discuss this further in Section 

The intermediate stellar mass galaxies (M» ~ 1O 9 I5 M 
to M» ~ 10 11 Mq) show a decreased rate of formation com- 
pared to the lower stellar mass galaxies. This manifests itself 
most clearly in the lowest redshift bin as a dip in the galaxy 
stellar mass function in the intermediate stellar mass range. 
This feature is also possibly present between redshifts of 
z = 2 — 3, although less obvious. We discuss this feature in 
Section \5\ 

As noted before we find a general trend for a is to 
increase at higher redshift, as shown by the black line in 
the middle panel of Figure [2] We also find that our values 
of a are more negative (therefore steeper) than found in 
previous studies. The middel pane l of Figure [2 ] shows that 
our results for a are steepe r than Eisn er et al ] (|2008h and 
iPerez-Gonzalez et al. I (l2008h. This is also the case for other 
studies, including IFontana et al. I (|2004l) (held a constant at 
-1.27/ - 1.36 from z = 1 - 2) and IFontana eta] (|2006h 
(a = -1.27 to -1.47 from z = 1.15 - 3.5) . We do however 
find a similarity between our results and iKaiisawal (|2009| ) 
who explore a similar redshift range and stellar mass depth 
to us. We find that <f>* is decreasing at higher redshift, as 
demonstrated by the straight line fit to the right panel of 
Figure [2] The overall decrease in cj>* represents the overall 
decrease in number density, as is expected at higher redshifts 
since fewer galaxies have had time to form. 



3.2 Blue/Red Mass Functions 

To obtain the galaxy stellar mass functions for the blue and 
red galaxy populations we divided the galaxies by colour. 
This was done by dividing the sample using the red sequence 
equation, 

(U - B) = -0.032(M S + 21.52) + 1.284 - 0.25, (4) 

from lWillmer et al.l (120061 ). modified for the AB magnitude 
system. The equation applies at redshift of z — 1.0, hence we 
modified Equation [4] to account for redshift evolution. This 
was done by using the redshift evolution of the luminosities 
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Figure 3. The division of the red and blue galaxy populations. The solid black line represents the colour cut used to divide the galaxies 
into their respective populations corrected for redshift evolution. 



Redshift Range M* 0(xlO~ 4 ) a 



1.0 to 1.5 


Blue 


11.51 ± 0.15 


3.07 ± 1.35 


-1.48 ± 0.06 




Red 


11.41 ± 0.16 


2.88 ± 1.09 


-0.60 ± 0.19 


1.5 to 2.0 


Blue 


11.06 ± 0.14 


6.94 ± 3.23 


-1.27 ± 0.13 




Red 


11.13 ± 0.12 


8.42 ± 1.69 


0.03 ± 0.30 


2.0 to 2.5 


Blue 


10.68 ± 0.21 


14.19± 10.16 


-1.38 ± 0.26 




Red 


10.58 ± 0.24 


4.79 ± 1.75 


0.82 ± 0.63 


2.5 to 3.0 


Blue 


10.67 ± 0.24 


10.21± 9.40 


-1.60 ± 0.29 




Red 


11.11 ± 0.29 


1.92 ± 0.85 


-0.06 ± 0.54 



Table 3. The values of the parameters from the Schechter fit to the blue and red galaxy population. 



and colours of galaxies from Ivan Dokkum fc Franxl l|200ll ). 
This allowed us to obtain a change in Mb and (U — B) 
at higher redshifts due to passive evolution. The change is 
then included in the equation so that the right hand side of 
Equation [4] becomes 

- 0.032(M S - AMs + 21.52) + 1.284 - 0.25 + A(U - B).(5) 

We then apply the cut so that if (U — B) is greater than 
Equation [S] the galaxy is red, and if (U — B) is less than 
Equation[5]the galaxy is blue. The results can be seen in Fig- 
ure(3] where the red/blue triangles/squares are the red/blue 
galaxies and the black line is the colour cut. 

After applying this colour cut we split the blue and red 
galaxy populations into the same redshift bins as before, 
and the number densities computed in the same way for each 
colour. Also the errors on these mass functions are computed 
exactly as before using our Monte Carlo approach (Section 



13 . 1 . 3 p . The resulting colour stellar mass functions are shown 
in Figure [4] Here the blue and red solid lines represent the 
Schechter fits to the blue and red galaxy populations, and 
t he solid black curv e is the local galaxy stellar mass function 
of lCole et ail (|200lh . The green points are the points that are 
not included in the fitting of the blue population. Both the 
blue and red stellar mass functions are fit using equation [3] 
In this case we see no evidence for a constant M* , and thus 
we leave all parameters free when fitting. This is discussed 
further in Section [5] in terms of implied evolution. 



3.2.1 Inspection of the Blue/Red Mass Functions 

The colour stellar mass functions in Figure [4] like the to- 
tal stellar mass functions, show stellar mass downsizing. We 
see that even in the redshift range z = 2.5 — 3, red and 
blue galaxies at M» > 1O 11,5 M0 are present with number 
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Figure 4. The colour divided galaxy stellar mass functions. The red and blue points represent stellar mass functions of the red and blue 
galaxy populations as defined in Section 13,21 The green points are the points not included in the Schechter fit . The solid red and blue 
curves are the Schechter fits. The solid black curve is the local galaxy stellar mass function of lCole et al. | fcOOfT ). The solid pink vertical 
line is the theoretical mass limit of this survey (see Section lij.l.lH . 




Figure 5. The Schechter fit parameters of the blue and red galaxy stellar mass functions. Left: The results of the fitting for M*. The 
blue/red points are M* for the blue/red galaxy populations when all of the parameters are free. Middle: The result of the fitting for a. 
The blue/red points are a for the blue/red galaxy populations when M* is held constant. Right: The result of the fitting for 0*. The 
blue/red points are a for the blue/red galaxy populations when M* is held constant. All the red points are offset by 0.05 in redshift for 
clarity. 
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densities very close to the local value. In the redshift bin 
z = 1 — 1.5 we see that both the red and the blue galaxies 
with Af, > 10 U S M Q are nearly fully in place, whereas the 
low stellar mass galaxies have not formed as quickly. 

The low stellar mass end shown in Figure|4]is dominated 
by blue galaxies at z < 3, and their number densities are 
close to the local stellar mass function out to a redshift of 
z = 3. This means that the steepness of the slope in the total 
stellar mass function (Figure [TJ is dominated by low mass 
blue galaxies. For the blue galaxies themselves a remains 
roughly constant and <f>* is constant within the error bars. 
This is shown in Figure The parameter Af* on the other 
hand, shows a general decline unlike what we see in the total 
stellar mass functions. 

For the red population, the fitted value of a shows a 
general increase (less steep) at higher redshifts, which is also 
the case for the parameter Af*. Contrary to this, the param- 
eter <f>* shows very little variation. Unfortunately, we do not 
have good number statistics for the stellar mass functions of 
the red population (in the range z — 2.5 — 3 there are 26 red 
galaxies compared to 639 blue) . We also have large errors on 
the high mass galaxies over z — 2 — 3, and hence the fit to 
the red stellar mass functions do not provide robust results. 
To this end, we cannot make any strong conclusions about 
the evolution of the red galaxy population. 

3.3 The Mass Functions of Star Forming and Non 
Star Forming Galaxies 

Having previously examined the colours of the GNS galaxies 
we next investigate the differences or similarities between 
colour and star formation selected stellar mass functions. 
We use the star formation rates calculated as described in 
Section ^. 4l for the sample between z = 1.5 — 3 and for Af* > 
10 9,5 Mq. To divide the galaxies into passive and evolving 
populations we use the star formation rate divided by stellar 
mass of the galaxy to calculate the time it would take for 
a galaxy to double in size, we call this tdoubie- Using the 
Hubble time (th) at the redshift of each galaxy we obtain a 
measure of how fast a galaxy is forming by calculating tf orm , 
where 



^form ^double/^/i • 



(6) 



We averaged over all the values of tf or m to cut our sample 
into two distinct populations. Those above an average of 
(tform)=0.1 are considered non-star forming, those below are 
considered star forming. The resulting stellar mass functions 
can be seen in Figure [U 

The general trend in Figure [6] shows an increase of 
both populations at high stellar mass. In this region Af* > 
lO n Af we find the non-star forming galaxies dominate. At 
the low stellar mass end we see that the star forming popula- 
tion dominates over the non-star forming galaxies. The slope 
of the low stellar mass end is very steep with this steepness 
decreasing over time, even within the large error bars. 

3.3.1 Comparison Between the Blue/Red and the Star 
Forming/Non Star Forming Mass Functions 

Over plotted on Figure [5] are the Schechter fits for the blue 
and red galaxy populations, shown as the blue and red dot- 
ted lines. We find that over all the redshift ranges the star 



forming/non-star forming galaxies match the blue/red pop- 
ulations within the errors plotted. However there is some 
slight disagreement at the high mass end of the mass func- 
tions in the range 2 = 2 — 3. The highest mass star forming 
galaxies are better represented by the red Schechter func- 
tion, this is likely the result o f the presence of dusty s tar 
formation (Bauer et al. (2011). iGriitzbauch et al] (|201lh . 



4 INTEGRATED STELLAR MASS DENSITIES 
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Figure 7. The integrated stellar mass density calculated for 
each redshift bin using the integration of the Schechter func- 
tion. The black c ircles show the r esults of this work. The pink 
squares are from iDickinson et al.l j2003T) . the dark blue crosses 
are f romlPe rez-Go nzalez et al,l 1 1200811 , the light blue triangles are 
from lElsner et al l [|2008l l and the green and purple diamonds are 
from lDrorv et all ll2005t) .The black dashed dot line is a prediction 
of th e stellar mass densi ties from the integrated star formation 
rate jWilkins et all [20081. The final redshift point is plotted as 
an open circle as it is not considered in the analysis as discussed 
in Section l3.ll 



We integrate over the Schechter function to get the in- 
tegrated stellar mass density between z — 1 — 3. This is the 
total amount of stellar mass contained within galaxies, in 
a given redshift range, per comoving volume. The results 
of this calculation are shown in Figure [7] To calculate this 
we perform a numerical integration between the limits of 
Af» = 10 12 Af Q and Af, = 10 7 Af Q . We extended the integra- 
tion beyond the lower mass limit of our survey by extrapo- 
lating the Schechter function to masses beyond those which 
we fit. 

Our stellar mass densities are generally higher than 
what has been found in previous work, as shown in Fig- 
ure [7] due to our steeper values of a. The black dashed dot 
line shows the stellar mass density history obtained from 
the integration of the instantaneous cosmic star formation 
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Figure 6. The galaxy stellar mass functions of star forming and non-star forming populations, defined by Equation Jg}. The pink 
diamonds show the non-star forming galaxy populatio n and the blue half circles are star forming galaxies. The solid black curve is the 
local galaxy stellar mass function of lCole et all <200ll) . The solid pink vertical line is the mass limit of this survey. The red and blue 
dashed lines are the Schechter fits for the red and blue galaxy populations plotted for comparison (see also Figure l4l. 
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Table 4. The values of the parameters from the Schechter fit. M* is the mean value from fitting with all parameters free. <f>* and a are 
the result of the Schechter fit with M* held constant. 



history computed by IWilkins et ail 1.2008), using the same 
IMF as in our work. They show that at z > 0.7 the line does 
not agree well with previous work which inferred from the 
integration of the stellar mass at these redshifts. The stel- 
lar mass densities that we calculate are higher than what 
has been found in previous works, but are still slightly sys- 
tematically lower than the integrated star formation history 
line. 

We also computed the stellar mass density with a bright 
lower limit of M» = 1O 9 ' 5 M0 and a faint lower limit of 
M t = 10 3 Mq. The bright limit allows us to compare our 
computed stellar mass densities with the mass we can actu- 
ally observe in all mass bins, and the faint mass limit allows a 
better comparison with the integrated cosmic star formation 
history. We find that when we adopt these bright and faint 
limits we see very little difference between these the com- 
puted stellar mass densities. Even in the range z — 2.5 — 3, 
where the largest difference is found, the stellar mass den- 
sity from both the bright and faint limit still lie within the 
errors of the original stellar mass density. 

The stellar mass densities found in this work can be 



compared to the local value of ICole et alj (|200lf ). We find 
therefore that by z ~ 1.25, 40lg°% of the stellar mass in 
the local universe has formed and by z ~ 1.75, 4lj^J % has 
formed. The growth of stellar mass does not alter dramati- 
cally after z — 2. At z — 2.25, 30lg°% has formed, and we 
see that at z = 2.75, 30±8% is in place. The final stellar 
mass density at z = 3.25 is 14 ±3% of the local value, al- 
though the results for this redshift range are not as robust as 
the others, as discussed in Section \3. II Thus, for 1 < z < 3, 
we see that the stellar mass density changes only slightly, 
hence there is very little change in stellar mass over this 
time period. We also note the stellar mass forms quickly as 
roughly one third of the stellar mass is already formed when 
the Universe is only 2.3 Gyrs old. 

Figure [8] shows the stellar mass density for the blue and 
red galaxy populations, with the total stellar mass density 
shown in black. We see a large decrease in the mass den- 
sity for the red galaxies, yet for the blue galaxies the stellar 
mass density is consistent with being constant over this red- 
shift r a nge. This is see n in ot her st udies such a s Borch et al.l 
( 2006). lArnouts et al.l (|2007T ) and lllbert et"aH l|2009bl ) who 




Figure 9. The ratio between the local galaxy stellar mass functions of lCole et al. and the number densities of the total galaxies in 

this work. The green line is the local stellar mass function over the local stellar mass function. The dotted line shows the local stellar mass 
function with varying <f>* and the dashed line shows the local stellar mass function with varying a. The final redshift bin is represented 
by open circles as it is not included in the bulk of our analysis as discussed in Section 13.11 



combined find very little evolution in the blue stellar mass 
functions as far out as redshift of z ~ 2. For the red popu- 
lation the numbers themselves are not as robust due to the 
poor fitting, but the different evolution of the red galaxies is 
clear. This growth in stellar mass density in the red popu- 
lation is co nsistent with other works such as lArnouts et al] 
l|2007D and lllbert et~al] (|2009bT ). 



5 DISCUSSION 

5.1 Comparison With Other Surveys and Data 

There is a large amount of research exploring the stellar mass 
function over a large range in redshifts, from z = — 6. What 
is unique about the GNS, and this work in particular, is the 
depth and redshift range provided by the data. A large num- 
ber of surveys (e.g. COSMOS, zCOSMOS and COMBO-17) 
previously examine the stellar mass function at a lower red- 
shift range than the GNS. For z < 1 these surv eys probe a 
large stellar mass range e.g. iPozzetti et all |2009) , with mass 
limi ts of M, ~ 1 8 - 5 M o at z = to 1O 1O M at z = 0.75 
and lDrorv et al.l (|2009T ). with mass limits of Af» ~ 10 8 - 5 M© 
at z = 0.2 to 10 9 M© at z = 1. 

For z — 1 — 1.5 our stellar mass functions are com- 
plete down to M* ~ 1O 8 5 M , and for z = 3 - 3.5 we 
find we are complete to M* ~ 1O 9 ' 5 M0. Other work, which 
probes a similar redshift range, does not probe the faint 
end of the mass function as low as we do here. For exam- 



ple iPerez-Gonzalez et al] l|2008l) has a mass limit of M* ~ 
10 9 - 5 M W at z = 1 - 1.3 and the mass limits of iKaiisawal 
(|2009h range from M* ~ 10 9 M Q to 1O 1O M over z = 1-3.5. 
We are thus able to examine stellar mass evolution at z > 1 
as a function of mass for the first time. 



5.2 The High Mass Galaxies and Downsizing 

In all of the stellar mass functions presented in this work, 
we find that the high mass galaxies have formed earlier than 
the low mass galaxies. For the total galaxy stellar mass func- 
tions, as shown in Figure [T] galaxies whose stellar mass 
is M» > 10 11,5 Mq are close to the local value as early as 
z = 3 — 3.5 (for M* ~ 10 11,5 Mq the number densities are 
3lli4 % of the local value at this redshift). These galax- 
ies have reached the local number density by z = 1 — 1.5 
(94^3 9 % the local value). Downsizing conventionally says 
that the high mass galaxies have stopped star forming be- 
fore the low mass galaxies. However what we see here is 
galaxy "mass downsizing", where by the mass of the high 
mass galaxy is formed before that in the low mass galaxies. 
This is occurring at z > 3 unlike star formation downsiz- 
ing which starts at z < 1.5 (Bauer et al. 2010). This means 
the majority of the stellar mass of a galaxy is already in 
place before the star formation largely stops. This suggests 
that it is not the star formation making the galaxy mas- 
sive, instead it is some other process such as mergers, which 
we do see happening often at high redshift l|Conselice et al] 
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Figure 10. The total Schechter functions (green dashed line). Over plotted is a double Schechter function (black dash dotted line) 
calculated using the blue and red fit added together. 



120031 . IConselice ct al. 2006 and Conselice et al. 2008). This 
mass downsizing effe ct is also seen for the halos of galaxies 
l|Foucaud et a]||2010l V 

Figure [5] also shows the rela tive ratio between the local 
galaxy stellar mass function of ICole et al.l (|200lr ) at z = 
and the number densities calculated in this work. The 
green horizonal line is how the ratio would appear if all the 
number densities were at the local value. The black line is 
a linear fit to the points from the lowest to highest ratio, 
hence its steepness correlates to how dominant downsizing is 
in that redshift range. We find that although at z = 2 — 2.5 
the gradient of the line is very steep, there is a generally 
decreasing trend with lower redshift. 

We find similar results for the red and blue mass func- 
tions in Figure[4] with the red population of Af* ~ lO n ' 5 M0 
galaxies being very close to their local number density in the 
range z — 2.5 — 3 (15lg 8 % for galaxies in this redshift range). 
The evolution of the high mass blue galaxies towards the lo- 
cal value is slower than that of the red. We find that the 
blue galaxies with stellar mass M 4 > 1O 11,5 M0 are mostly 
present by z = 1 - 1.5 (56tll% for M* > 10 n ' 5 M Q for 
galaxies with in this redshift range). This is because these 
blue galaxies are likely to evolve into red galaxies at some 
stage. We also find downsizing in the star forming and non- 
star forming mass function, where the star forming/ non-star 
forming populations behave very similarly to the blue/red 
populations. 



5.3 The Intermediate Mass Dip 

Figure [9] shows various lines indicating the ratio between 
the local stellar mass function and the local stellar mass 
function with slightly varying parameters (the coloured and 
dashed/dotted lines). For each line only one parameter is 
changed at one time. We see that between z — 1.0 to 1.5 
galaxies in the mass range M* ~ 10 10,5 M Q - 1O 9,5 M , have 
generally higher ratios, hence they have generally lower num- 
ber densities relative to galaxies with M* > 1Q 11 Mq and 
M» < 1O 9 5 M . This dip feature also seems to be present 
for z > 2, but is much weaker and shifted to the stellar mass 
range M» ~ 10 115 M Q - 1O 1O - 5 M . The various coloured 
lines show the features of the mass function is not just an 
effect of the changing of the shape of the Schechter function. 

The dip in the total stellar mass function has been seen 
before and can be explained by the differential evolution 
of the blue and red galaxy populat i ons. T his is consistent 
with studies such as IPozzetti et all (2009 and lllbert et al.l 
(|2009bh who find a dip in the total stellar mass function that 
can be fit by the combination of these two population. We 
also find this to be the case here, as shown in Figure [TO] We 
find that a double Schechter function (the black dash dotted 
line), computed using the parameters of the blue and red 
fits, matches the form of the total Schechter function for all 
galaxies. The only large discrepancy occurs at the high mass 
end in the range 2 < z < 2.5, and this corresponds with a 
poor fit to the red mass function. 

iDrorv et al.l {2009) find that out to a redshift of z = 
1, a dip in the stellar mass function is observed not only 
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Figure 8. The integrated stellar mass density calculated for each 
redshift bin using the integration of the Schechter function for 
the red and blue galaxy populations. The black circles are the 
total stellar mass density plotted for reference. The black points 
are offset in redshift by -0.05 and the red points are offset in 
redshift by +0.05. The final redshift point for the total population 
is plotted as an open circle as it is not considered in the analysis 
as discussed in Section 13. II 



in the total but also in the red and the blue stellar mass 
functions for M* ~ lO lo A/0 galaxies. They find that the 
blue stellar mass functions show a dip feature that appears 
to become more pronounced between z = 0.2 and z = 1. The 
also observe a dip in the red mass function that shows the 
opposite trend, becoming more pronounced at lower redshift. 
The fact that we do not see this feature in either the blue or 
the red mass functions tentatively suggests that any causes 
of the dip are perhaps not present beyond 2 = 1 for these 
colour selected populations. 



5.4 The Low Mass Slope and M* 

Figure [2] shows a general steepening of a from low to 
high redshift. It is one possibility that the low stellar 
mass galaxies are formed early on, then undergo mergers 
where we k now that these merger events happen at an 
early epoch (IConselice et al.ll2005l . lConselice et al. I l200d and 
IConselice et all 120081 ). These mergers form higher stellar 
mass galaxies and hence the the number densities of lower 
mass galaxies would decrease over time, causing a to become 
less steep, as is observed. They could also fade, creating a 
less steep mass function. 

In this work we generally find a steeper a than pre- 
vious investigations of the stellar mass function at high 
reds hift. For example Figu r e [2] s hows we are steep e r tha n 



both|P erez-Gonzalez et all (|2008l 'l an d lElsner et al 



We also find a steeper slope than iFontana et al 
(a = —1.36 over z = 1 — 2) and I Font ana et al 



2008) 



2004) 



2006) 



(a ~ —1.3 at z ~ 1.15 then steepens toa~ —1.5 at z ~ 3.5). 
In this case, the higher stellar mass densities discussed in 
Section U are most likely the result of the steeper slope. The 
depth of the GNS allows us to probe deeper into the low 
stellar mass region, allowing us to uncover the steeper slope 
and higher stellar mass density. We test this by recalculating 
stellar mass densities with less steep values for the param- 
eter a, but keeping the remaining parameters the same as 
found in this work. We find that by increasing a by 0.1 in 
the redshift range z = 2.0 — 2.5 there is a ~ 16% decrease in 
the integrated stellar mass density. This decrease becomes 
~ 27% when alpha is increased by 0.2 in the same redshift 
range. IPerez-Gonzalez et alj (2008) find a = —1.26 in the 
redshift range z — 2.0 — 2.5 and hence compute a stellar 
mass density of 7.87±0.09. Using this value of a we compute 
a stellar mass densit y of 8.06 ± 0.19, which is w ithin error 
of the value found by I Perez- Gonzalez et al.l |2008l ) . Since we 
can reproduce the lower stellar mass densities with a shal- 
lower value of a, our higher stellar mass densities are likely 
a result of our steeper a, and additional mass we see at the 
low mass end of the stellar mass function. We find that the 
steepness, and hence the additional stellar mass we see, is 
dominated by the low stellar mass blue galaxies we find. In 
Section 13.1.31 we test the effect of the errors on our number 
densities and find we can still reconstruct the Schechter pa- 
rameters well, thus errors calculated in this work are ruled 
out as being the cause of our steeper a. 

For the blue mass functions in Figure [4] a is consistent 
with being constant across the whole redshift range. There 
is also very little evolution of the parameter <j>* , and this 
leads to the relatively unchanging blue stellar mass density 
seen in FigureJS] This is consistent with lllbert et al l (|2009bh 
who show that between z = 0.2 and 2 both the stellar mass 
function and the stellar mass densit y show very litt l e evo - 
lution. This is also in agreement with lVergani et al. I (|2008l ) 
who shows the parameters for the blue galaxies are virtually 
unchanged between z = 0.5 and 1.3. This suggests there is 
some process by which this stellar mass function is being 
replenished. 

We find that for the total mass function, the parameter 
M* stays roughly constant over time. This implies that the 
process which is increasing the numbers of galaxies at the 
high mass end, is doing so over the whole range of high 
masses above and around M* . However this is not seen with 
the blue galaxies, thus this effect is not part of the evolution 
of this population. 

There will be some degeneracy between the Schechter 
function parameters, and hence by holding M* constant we 
may be affecting the results of a and <j>. To test the validity 
of the assumption of a constant M* we use a Monte Carlo 
simulation to see how the parameters a and <f> vary with M* . 
To do this we took the values of M* we obtained from the 
first fitting described in Section [3. 1.21 and then varied these 
values between the extremes of the error on M* . This gave 
us a range of M* and we then repeated the fitting for each 
value in the range keeping a and <j> free, thus showing us how 
the other parameters vary. We also recalculated the stellar 
mass densities using Monte Carlo Schechter parameters, the 
parameters and stellar mass densities are also plotted in Fig- 
ure 1111 We only apply this analysis to the range z = 1 — 3 
as the final redshift bin is mostly excluded from the analysis 
as described in Section [3Tl 
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We find that although variations in M* does produce a 
spread in the parameter results we still see the same general 
trends. We still produce steep values for a, which in turn 
produce stellar mass densities that are higher than found in 
previous studies. We plotted histograms of the parameters 
and stellar mass densities and found that they do peak at 
the values found from fitting with constant M*. 



6 SUMMARY 

We construct galaxy stellar mass functions using 8298 galax- 
ies detected in the GOODS NICMOS Survey, a HST H 160 - 
band imaging survey centred around the most massive 
galaxies in the redshift range z — 1.7 to 2.9. We use the 
stellar masses, redshifts, colours and star formation rates 
calculated from the H-band imaging combined with the ACS 
BViz data to calculate stellar mass functions over this red- 
shift range. We calculate the mass limit of the survey and 
find that we are probing galaxy stellar masses down to 
M* ~ 1O 8,5 M0 in the lowest redshift bin, hence we are prob- 
ing an unexplored region of the stellar mass function whilst 
maintaining statistically meaningful results (~300 galaxies 
between z = 3 — 3.5). We correct for high mass bias, consider 
both redshift and stellar mass errors, and fit Schechter func- 
tions to examine the evolution of the parameters. We also 
measure the galaxy stellar mass functions for a blue/red 
selection and one based on star forming/non-star forming 
galaxy populations. We also calculate the corresponding to- 
tal stellar mass densities by integrating over the stellar mass 
functions at each redshift range. 

The major results of this paper are as follows: 

• We observe stellar mass downsizing in all of our stellar 
mass functions. For the total galaxy stellar mass functions 
we see that galaxies at M* > 10 11 Mq have reached nearly 
the local value by z — 3. Even in the highest redshift bin it is 
clear that the most massive galaxies are very close to being 
in place. The same is true for the red and blue galaxies. 

• At all redshifts the blue galaxies dominate the low mass 
end of the stellar mass functions. Between 2 = 1 — 3 the 
lowest stellar mass blue galaxies are close to the local value. 

• The total stellar mass function shows a dip feature in 
the intermediate mass range, which can be explained by the 
differential red and blue populations evolution. We see no 
such dip in the blue and the red stellar mass functions unlike 
some previous studies who have found this below z ~ 1. We 
suggest these features are not present at high redshift. 

• We find a generally steeper low mass slope for the total 
stellar mass functions than previous work due to the low 
stellar mass blue population that is probed by our deep data 
set. This results in a generally higher stellar mass density. 

• We find a generally higher stellar mass density, due to 
probing deeper into the lower stellar masses. 

By constructing the galaxy stellar mass function we are 
investigating the mass differentiation in galaxy formation, 
and the way in which mass drives evolution. More work must 
be done to fully understand the similarities and differences 
between the formation of the high mass galaxies that form 
first, the intermediate mass galaxies who seem to show a 
decreased rate of formation, and low mass galaxies who are 
dominated by the blue population. Large surveys such as 



CANDELS and future spectroscopic surveys that could be 
performed with telescopes such as JWST and E-ELT, will 
provide the quality of data required to continue advancing 
in this field. 
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Figure 11. The results of the Monte Carlo simulation to test the validity of the assumption of a constant M* . The black circles are the 
parameters a and 0, and the stellar mass densities found from the range of M* . The pink circles show the original a and and stellar 
mass density from the fitting described in Section l3.1.2l The right hand panel also includes results from literature. 
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